Preparation of NOON State Induced by Macroscopic Quantum Tunneling in an Ising 

Chain 
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In this brief report, we propose a possible way, theoretically and experimentally, to generate a 
NOON state of the two degenerate ferromagnetic ground states of the Transverse Ising Model. In our 
scheme we employ the macroscopic quantum tunneling (MQT) effect between the two degenerate 
ferromagnetic ground states to realize the NOON state. Our calculation about the MQT process 
is based on a higher-order degenerate perturbation method. After doing a transformation, the 
MQT process could also be treated as the hopping of individual virtual fermions in the spin chain, 
which will leads to an analytical description of tunneling process. The experimental feasibility for 
generating the NOON state is discussed in the setup of linear ion trap. 



The NOON states are a kind of special quantum states 
with two orthogonal component states in maximal super- 
position, where one is all spin-up and the other is all spin- 
down, such as \NOON) = (Iff • • • ft) ± III ■ • • W» /V2. 
The recent interests in production of NOON states 
mainly arise from the crucial role they play in quan- 
tum optical lithography [l| and quantum metrology [^-Q 
as well as quantum information processing Q, including 
quantum gating Q and precision measurement of phases 
0-Q . Many efforts have been devoted to preparation of 
NOON states theoretically and experimentally by using 
trapped ions [Iol - [l^ | , cold atoms [H, HH , atomic ensem- 
bles superconducting quantum circuits [IH, cavity 
quantum electrodynamics |17ll. and linear or nonlinear 
optical elements [j, 0, [H|-|2(|7 

In this Brief Report, we propose an idea to create the 
NOON states formed by two degenerate ferromagnetic 
ground states of a transverse Ising model by the mecha- 
nism of macroscopic quantum tunneling (MQT). In our 
case, the transverse Ising model is given by 



H = H + H! 



(1) 



where Hq 



^E^^k and H, = h'^af. In 



general, the transverse Ising model includes two quan- 
tum limits at absolute zero temperature: the ferromag- 
netic (FM) limit (J > h x ) and the paramagnetic limit 
(J < h x ), separated by the strength of the transverse 
field J c [2l[. In the present model, Hj is the pertur- 
bation term to induce quantum tunneling, as discussed 
later. The global phase diagram of the transverse Ising 
model at zero temperature, as shown in Fig. 1, has been 
well known [HlH, where quantum phase transition oc- 
curs from the ferromagnetic state with J > h x to the 
gapped spin-polarized state with J < h x at the quan- 
tum transition point J/h x = 1. As a result, for large J 



( J > h x ), the ground state owns a spontaneous non-zero 
magnetization in z direction due to its FM order. For 
small J (J < h x ), the ground state is the uncorrected 
spin state, also called the spin-polarized state. 

The key idea in our proposal is the quantum tunneling 
between degenerate ground states induced by a trans- 
verse perturbation field. In a macroscopic system, such a 
quantum tunneling is usually referred to MQT, which 
has appeared in a wide range of research fields, such 
as quantum oscillations between two degenerate wells of 
NH 3 , quantum coherence in one-dimensional charge den- 
sity waves, the ferromagnetic single domain magnets and 
biased Josephson junctions [23j]. In general, the MQT can 
be found in systems with two or more separated 'classical' 
states that are macroscopically distinct. For instance, the 
MQT in macroscopic systems at low temperature, such as 
the superconducting quantum interference device, might 
be feasible for observation [24j . 

There have been some useful theoretical methods 
to deal with the MQT, such as semiclassical WKB 
(Wentzel, Kramers and Brillouin) method [25[ and In- 
stanton method in path integration [26j | . In the present 
work, we study the MQT in a transverse Ising chain by 
the perturbation method [27l - [29| , based on the fact that 
the transverse Ising model in a finite size owns the degen- 
eracy on the ground state, which can be removed by the 
MQT process. This mechanism can result in production 
of NOON states in the transverse Ising model. 

To classify the degeneracy of the two ground states 
of H, we introduce two string operators Wx = Yl&f 

and W z = n°f E3, HI, which commute with fl, 



H ,W X 



H ,W Z 



o, i.e., 

= and are anti-commutated 
{W z ,Wx} = 0. Due to these 
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with each other, i.e. 
properties, we may represent Wx and Wz by two- 
component pseudo-spin operators t x and t z , respec- 
tively, as Wx — > t x and Wz ~> t z , for which we have 
r z III • • • II) = - III • • • II) ,r z Iff • • • ft) = |tt • • • Tt>, 
where |tt " " ' tt) an d III ■ ■ ■ II) ar( 3 two degenerate fer- 
romagnetic many-body ground states, as show in Fig. [T] 
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FIG. 1: The global phase diagram of transverse Ising model at 
zero temperature, where the quantum phase transition hap- 
pens as J crosses the critical value J c . Our model stays in 
the ferromagnetic limit ( J 3> h x ) of the phase diagram. 

In order to illustrate the MQT process in the Ising 
chain, we transform the Ising spin chain into the system 
of free fermions by a perturbative method [28]. In the ab- 
sence of Hi, the fermions in this model have flat bands 
with the energy spectrum of E(k) = 2J, which implies 
that there is a mass gap 2 J of the fermions in the trans- 
formed system, and the fermion cannot move in such a 
case. To change this situation, we may introduce the 
transverse perturbation Hj , which drives the fermions to 
hop in the model. 

In the perturbative approach, the perturbative oper- 
ators af can be represented by hopping terms of the 
fermions, af -> (ipjipi+i+iplipl^+h.c), where tpj (tpi) is 
the creation (annihilation) operator of a spinless fermion 
at the site i. In addition, we introduce a constraint for 
single occupation per site (i.e., hard-core constraint), 

(V>J) 2 |*) = I*) or $\*)=r/, i \V), 

where |\&) denotes a many-body quantum state. 

Therefore, by the perturbative method, the Hamilto- 
nian in Eq. (1) can be rewritten by spinless fermions 
as, 

H = 2J^2 $1>i - h X E + $W + h '°)- ( 2 ) 

Under the Fourier transformation, we may assume peri- 
odic boundary condition for the Ising chain, and obtain 
the dispersion of the fermion system in Eq. ^) as 

E(k) = V(2/i x cosfc + 2J) 2 + A{h x sinfc) 2 , (3) 

where k is the wave vector, and for simplicity we set 
the separation of two nearest neighbor spins to be unity. 
So the energy gap of the fermions is obtained asmf = 

2 J^J 1 — ij-, and thereby we may manipulate the disper- 
sion of the fermions by tuning the external held h x . In 



FIG. 2: (Color online) The hopping process of the fermions 
from the site i to 1. At site i, a pair of fermions are gener- 
ated by the operation of. One fermion (at right cross in the 
lower string) hops when one does an operation <rf +1 . Then the 
hopping process of fermions is driven by the transverse per- 
turbation field h x . We represent the fermions by red crosses 
on the string. 

the language of the free fermion system, during the tun- 
neling process, a pair of virtual fermions are created, one 
of which moves around the chain until annihilation by 
meeting another fermion [27M29I I . 

In the thermodynamic limit, the two ferromagnetic 
ground states of Eq. (JTJ are degenerate. However, there 
is very tiny energy splitting AE between the two ground 
states in a finite size chain. Since AE is much smaller 
than the energy gap in the fermionic excitation, we 
may safely ignore those excited states and only concen- 
trate on the MQT between the two degenerate ground 
states in our perturbative calculation. 

Now we employ a high-order degenerate perturbative 
app roach to calculate the MQT effect in our system [27i — 
|22j. Firstly we introduce the transformation operator 
Ui(0, —00) defined in units of h = 1, 

{7/(0, -00) =f exp(-i I H'j{t')dt'), (4) 

where H'j{t) — e lHat Hje~ tHot is for Hcisenberg picture, 
and T is the time-ordering operator. The transformation 
operator /7/(0, —00) in Eq. ((4]) can be further written as 

00 

17/(0, -co) |m) = ]T U { / } (0, -00) \m) , (5) 
3=0 

where 

?7j 0) (0,-oo) \m) = \m) , 

1>P 0) (0, -00) |m) = (— J-g-^y \m) (6) 
Eq — Hq 

( Iff . . . ff) \ 

and \m) = I j|| J. The element of the trans- 

formation matrix from the state |m) to |n) becomes 
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(n| Ui(0, -co) \m), corresponding to the energy, 

E= {n\ HUi{0,-oo) \m) = E + SE, (7) 

with E the eigenvalue of the Hamiltonian H for the 
cigcnstatc |m). 

When the MQT occurs between the two degenerate 
many-body ground states |ft ' ' • tt) an d 144- ' " 4-4-); wc 
may describe the process in spinless fermion language 
as that, a pair of spinless fcrmions are created, and then 
one of them travels around the chain and leads to a string 
operator. As a result, the dominant term in Eq. ([5]) is 
labeled by j = N — 1, where TV is the site number of 
the Ising chain. Considering the tunneling process, we 
obtain the perturbative energy as 

5E = {n\H I tfi(0,-oo)\m) (8) 

OO 

= (n\H I J2uP(0,-oo)\m) 

3=0 

= (n|iJr#r-%-oc)|m), 

where the operator Hi(j\ N ~ 1 \o, — oo) is proportional to 
a string operator Wx,z- 

To calculate the MQT between the two ground states, 
we consider a virtual fermion propagates around the 



chain, leading to the quantum state 



flipped as 



in • • • tt) 
in ■ • ■ id 



in • • • id \ _ ,( in • • • tt> A 



in • • • tt) 



in • • ■ id J 



(9) 



Using Eq. (jSJ), we may obtain the energy shift 8E as 



SE = U} 



(N) 



(tt ■ • • ttl Hj( 



1 



Eq — He 



H I ) N - i \U---U). (10) 



Due to the translation invariance of Eq. ([!} , we choose 
the site i as the starting point of the MQT process and 
wc have 



1 



En — Hi 



Hi) III • • • II) 



h x 



Eq — H 



°f)\U---U) 



En — Hn 



■)l*i>, 



(11) 



where is the excited state of the two fermions at 
the sites i and i — 1 with the energy En + 4J. Using 
H = (-Bo + 4 J) |*,), we have 

i—^-Hj) \u~- id = n(^-) m . 

En — Hq * J 



After the two fermions created, one of them moves 
along the Ising chain with following steps, 



h x 



En — Hn 

, h x 



E<) 2 i4-i •••id 



h x _ h x 



4 J Eq — Hq 



ux ux 

Eq-Hq + 

ux ux 

-4J n -4J n j/ ' 

where is the excited state with two fermions at sites 
i + 1 and i — By this way, the fermion moves around the 
chain step by step. When coming back to the starting 
point the fermion annihilates with the partner fermion, 
which implies that the system changes with the ground 
state flipped from ||| • • • ||) to |tt • • • tt)- So we finally 
get the energy splitting 

AN) 



AE = 2SE = 2U} 

= 2<tt---tt|#/( 



£fl — Ho 



H I ^- i \n---u) 



(hx\N ux 



(12) 



For the chain with open boundary condition, the tun- 
neling process is much different from our discussion above 
under the periodic boundary condition. In this case, the 
dominant tunneling of a single virtual fermion is from 
one side to the other with the tunneling splitting 



AE = 2SE 



x\N 



(2J) 



N-l ■ 



(13) 



By this high-order perturbative approach we have 
studied the tunneling splitting between the two degener- 
ate ground states. Fig. 3 shows the MQT calculated nu- 
merically by the exact diagonalization technique, which 
indicates that our theoretical results agree well with the 
numerical results. After adding the perturbation field h x 
to the Ising chain, we are able to generate two fermions 
and even drive one of them hopping around the chain by 
considering high-order perturbation terms. In this way 
we could manipulate the tunneling splitting of the two 
degenerate ground states. In addition, in the sense of 
'perturbative treatment', our approach can only be ap- 
plied to the cases under a small external field. Never- 
theless, the tunneling induced by the transverse distur- 
bance field h x in our model makes it possible to realize a 
string NOON state of the many-body state components 
|tt---tt) and W ■••44-)- 

Now we discuss how to generate a NOON state in a 
linear trapped- ion system. We may employ S±/2 (jnj = 



4 
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FIG. 3: (Color online) The energy splitting between the two 
degenerate ground states in an external field along ^-direction 
(J = 1). Here L denotes lattice the number. 

-1/2) and D 5/2 (rrij = -1/2) of 40 Ca+ [31] as the spin 
states |t) and \\.) in above transverse Ising model, respec- 
tively, where the lifetime of D 5 / 2 of 40 Ca + ion is longer 
than 1 sec. To carry out the present scheme in the ion 
trap, one should firstly confine and cool ./V 40 Ca + ions 
in a linear Paul trap, and make them located equidis- 
tantly in an array. Then we could generate a transverse 
Ising model by coupling the internal states of the ions 
and their collective motional phonon modes by choosing 
suitable values of the trap frequency and Rabi frequency 
of the system [l(| [IH, l32| . Experimentally, such a trans- 
verse Ising model has been simulated by trapped ions 
[33[ | . Furthermore, several other methods are reported in 
trapped- ion system to create the \NOON) states, such 
as the methods in [13, and more recently in [l2[ • 

In our case, after introducing the nearest-neighbor 
Ising interaction between neighboring ions, the whole 
Ising chain should be polarized to one of the two de- 
generate ferromagnetic ground states by irradiating the 
ion-chain with specially polarized lasers. For example, we 
can initialize the Ising chain to its ferromagnetic ground 
state |44 ' ' ' 44) which corresponds to the N spin-down 
of the chain. The following task is to realize the trans- 
formation operator [7/(0, —00) in the system. Note that 



the dynamical evolution caused by the small transverse 
perturbation field h x can be expressed as N — 1-th or- 
der perturbation of [7/(0,— 00), which is proportional 
to the action of e lAWx . As we had mentioned before, 
the purpose we introduce h x is to induced the tunnel- 
ing process between the two degenerate ferromagnetic 
limit ground states |tt " " " TT) an£ l 144 ' ' ' 44)- It implies 
that the transverse field h x is just a little disturbance 
compared to the big Ising coupling J along z-direction. 
Based upon the consideration above, we may globally flip 
the whole ions in the chain in an adiabatic way to gener- 
ate the perturbative string operator Wx = W erf , which 

i 

can be complemented by applying a Gaussian beam with 
large width experimentally. This key step takes time 

of T = ~KE = { 2{h*)" h & where AE can be calcu- 
lated using Eq. ([7]) by high order degenerate perturbation 
method. If we take a entanglement involving 14 ions as 
an example the time we need for the Gaussian beam 
is on the order of r ~ 10~ 12 s. The ferromagnetic domain 
wall in the chain implies that two fermions are created 
and every fermion has a 2J excitation energy which is 
much greater than AE. In this way, we have achieved a 
NOON state \NOON) = a |ft • • ■ tt> + ^ 144 " • 44) 
in the trapped ions under this transverse Ising model, 
the three parameters: a, (3 and (f> could be determined 
by experimental measurements. 

In summary, we have proposed a scheme to generate 
the NOON state of the two degenerate ferromagnetic 
ground states in the transverse Ising model through MQT 
process. This processes could be effectively manipulated 
by the perturbative field which is experimentally feasi- 
ble in ion-trap systems. We think that this Ising string 
NOON state may have potential application in quantum 
computation and even in topological quantum computa- 
tion if we could introduce topological protection to this 
Ising string. 
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